Lower limit on the achievable temperature in resonator-based sideband cooling 



00 
O 

o 



M. Grajcar, 1,2 S. Ashhab, 1,3 J.R. Johansson, 1 and Franco Nori 1,4,5 

frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi 351-0198, Japan 
2 Department of Experimental Physics, Comenius University, SK-84248 Bratislava, Slovakia. 
''Physics Department, The University of Michigan, Ann Arbor, 



USA 



Michigan 48109-1040, 
■* Center for Theoretical Physics, Physics Department, 
Applied Physics Program, Center for the Study of Complex Systems, 
The University of Michigan, Ann Arbor, Michigan 48109-1040, USA 
5 CREST, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 

(Dated: February 11, 2008) 

A resonator with eigenfrequency ui T can be effectively used as a cooler for another linear oscillator 
with a much smaller frequency uj m >C 0J r . A huge cooling effect, which could be used to cool 
a mechanical oscillator below the energy of quantum fluctuations, has been predicted by several 
authors. However, here we show that there is a lower limit T* on the achievable temperature, given 
by T* = TmCJm/iOr, that was not considered in previous work and can be higher than the quantum 
limit in realistic experimental realizations. We also point out that the decay rate of the resonator, 
which previous studies stress should be small, must be larger than the decay rate of the cooled 
oscillator for effective cooling. 
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PACS numbers: 85.85.+j, 45.80.+r 



INTRODUCTION 



Recently, a tremendous experimental effort has been 
devoted to the task of cooling mechanical oscilla- 
tors below the energy of quantum fluctuations. In 
spite of many experimental improvements, the quan- 
tum limit has not been achieved* 1 * 2 -' 3 -'^ Several papers 
that propose cooling mechanisms using electromagnetic 
(rf, microwave or light) resonators&£&2 or other cool- 
ing mechanism o 10 ! 11 ' 12 ! 13 ] 14 ! 15 to fulfill this task have ap- 
peared recently. These papers predict an enormous cool- 
ing effect. However, they do not explicitly state that 
there is a lower limit on the achievable temperature, as- 
sociated with the ratio between the frequencies of the 
coolant and cooled oscillators, that cannot be overcome 
and can play an important role for realistic experimen- 
tal realizations. Moreover, some formulas that appear 
in the literature can give temperatures below this limit, 
which will be described in more detail below. This lower 
temperature limit can be important for the most feasible 
designs using rf or microwave resonators. 

The electromagnetic resonators can be easily imple- 
mented on-chip beside a nano-mechanical oscillator and 
kept at low temperatures. Such structures, also known 
as MEMS and NEMS (micro and nanoelectromechani- 
cal systems), have already been realize d 16 ' 17 : 18 , achieving 
high frequencies in the GHz range (ui m ~ 10 GHz) but 
with small quality factors Q m ~ 10 — 500. Nevertheless, 
they can be used as sensitive elements for weak force- 
displacement detectio n 19 : 20 and they have been proposed 
as possible qubitsj 21 : 22 For such systems, the frequency of 
the basic mode of the NEMS becomes comparable to the 
resonance frequency of the electromagnetic resonator uj r . 
In this case, the temperature limit proportional to u> m /u> r 
which was negligible for optical-frequency coolers, can de- 
termine the lowest achievable temperature T*. The aim 



of this work is to call the attention of experimentalists 
to this fundamental limit which could help them design 
more effective cooling systems. 



II. 



SEMI-CLASSICAL APPROACH 



For the sake of simplicity, we will consider a RLC tank 
circuit (the results can be applied to any electromagnetic 
resonator, such as a transmission-line resonator, cavity, 
Fabry-Perot resonator, etc.). A mechanical oscillator is 
coupled to the capacitor such that the capacitance de- 
pends paramctrically on the displacement of the oscilla- 
tor. Such a system was thoroughly analyzed in Ref. |23| . 
and we only briefly introduce the equations of motion 
here. If the mechanical oscillator is a part of one of the 
capacitor electrodes, the capacitance C(x) » Co(l — x/d) 
depends on the displacement x of the oscillator from the 
equilibrium position, where Co = eS/do is the capaci- 
tance at x — and d is the renormalized distance be- 
tween the electrodes d — do/n. Here k is the coupling 
constant between the mechanical oscillator and the RLC 
circuit, and it can be expressed as the ratio between the 
mechanical oscillator capacitance C m , which depends on 
the oscillator displacement, and the total capacitance Co 
(we consider the case C m <C Co), i.e. k = C m /Co- If 
the RLC tank circuit is pumped by a microwave source 
Vp = Vpocosujpt, the voltage between the capacitor's 
electrodes is Vq = aj r V p o/T r , and the Coulomb energy 
of the capacitor depends on its capacitance which, in 
turn, depends on the oscillator displacement. Thus, the 
electromagnetic resonator and mechanical oscillator (i.e. 
cantilever) can be described by a system of differential 
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equations of two coupled damped linear oscillators 
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where T r ^ m are damping rates, uj r ,m are angular frequen- 
cies, Vf is a fluctuating voltage across the capacitor, Ff 
is a fluctuating force acting on the mechanical oscillator 
with mass M, and Q(t) = q p (t) + q/(t) is total charge 
charge on the capacitor. The Eq. ([2]) is nonlinear but 
can be linearized keeping in mind that we are interested 
to calculate charge fluctuations qj(t) which are much 
smaller than charge oscillations q p (t) driving by coher- 
ent microwave source. It is convinient to express qf{t) 
and Vf(t) in terms of quadrature amplitudes 

9/M = Qc(t) cosujpt + q s (t) smujpt , 
V/(t) = V c (t) cos ujpt + V s (t ) sin u p t , 



and rewrite Eqs. CQ),© in the dimensionless variables 



<7c, s 


= Qc,a/y/ C frhJ r , 


X 




T 




CO r ^p 


= LO r p/id m , 


f 


= fm/^m j 


f r 


= T r /2ui m , 


V c , s 


= u) r V CjS -\/Co/fkJ r , 


Ff 


= Ff/y/MultkOm , 


Vo 


= uj r Vo^/Co/4:MLo m Lj r d 2 



Here T r and T m are the temperatures of the electromag- 
netic resonator and mechanical oscillator, respectively. 
Using the slowly- varying-amplitude approximation^ and 
considering Langevin fluctuating forces caused by quan- 
tum noise^^ 



VcAt) = 



I LT r — -coth — — 
2 \2k B T r 



Ff(t) = y MT m tiuj m coth 
Eqs. ©,© read 

dq 



2k B T„ 



^L = -Aq(T)+F(r) 
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(4) 



F(r) 



V 



IV coth(huj r /2k B T r )^ c (r) 
f r coth(/iw I ./2A;B7r)£s(' r ) 

f m coth(7kj m /2fc B T m )£ m (Y) 




(5) 



q = (q c ,q s ,v,x), v = dx/dr and fj = uj p — Cb r . Here 
T r and T m are the base temperatures of the resonator 
and mechanical oscillator, respectively. Thus, we have 
a system of coupled Langevin equations^ which allow 
us to calculate the stationary covariance matrix defined 
as (7 = (qq T ) s for r — > oo. The diagonal terms of the 
covariance matrix determine the mean squared values of 
the vector components q. For example, = (v 2 ) s is 
the normalized mean squared velocity of the mechanical 
oscillator. The covariance matrix can be determined from 
the system of linear equations 

A a + aA T = B , 

where B is a correlation matrix defined as 
(FiirfFjiT 1 )) = Bij5(T -t'). If the fluctuating forces 
are uncorrelated, i.e. (£x(t)£ x ' ( T ')) = 5 x ,x'5(t — r') (here 
x, x' stand for c, s or m), B takes the form of a diagonal 
matrix with elements 



Bu — 



/ f r coth(huj r /2k B T r ) \ 
r r coth(huj r /2k B T r ) 
f m cotYi(ftw m /2k B T m ) 




(6) 



The mean value of energy of the mechanical oscillator 
fluctuations is 



f 



(7) 



Now, one can easily calculate the effective temperature 
of the mechanical oscillator from the definition relation 
for T* 



1 / t\JjJ m 
<Tvv = — cotn 

2 \2k B T* 



(8) 



As we will see later, the most appropriate parameters 
for cooling purposes are fj = —1, 2r m r r <C Vq <C 1. In 



this limit and for T r 
expressed as 



r m < r r 



a DD = - coth 



huj r 

2k B T r 



r r 



coth 



(9) 



Thus, the lowest temperature of the mechanical oscilla- 
tor is limited by the first term if the second term is made 
negligibly small by sideband cooling. As a matter of fact 
this term simply shows that even in our semi-classical 
approach we cannot 'cool' the mechanical oscillator be- 
low the zero-point energy which is consistent with the 
Heisenberg uncertainty principle. Indeed, it follows from 
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FIG. 1: The cooling factor T^/T m as a function of the normal- 
ized pumping amplitude Vo of the noiseless microwave source 
for u r = 10 3 , f m = 10~ 5 , f r = 10" 1 and k B T m = k B T r > 
huj r , m calculated numerically (circles) and from Eqs. (JSj,® 
(solid line). 



Eqs.O and © that the energy saved in the mechanical 
oscillator is 



■ coth 



fvjj r 
2k B T r 



(10) 



In this limit, the effective temperature of the mechan- 
ical oscillator takes the simple form 



UJ r 



(11) 



The cooling factor Tj^/T m as a function of the normalized 
pumping amplitude Vq is shown in Fig. [T] Even though 
this result was derived within semi-classical physics, the 
same limit can be obtained using the quantum approach, 
as we shall show below. 

Here we should emphasize that the temperature of the 
resonator T r is usually much higher than the ambient 
temperature if the resonator is heavily pumped by the 
microwave source. This is caused by the phase noise of 
the microwave source, which is directly proportional to 
the output power. Microwave sources are characterized 
by the single sideband noise spectral density^! 



L(Aw) = 10 log 



Si- 



lP 



(12) 



where U^ w 



is the mean-square voltage of the microwave 
source and Sy is the spectral density of the volt- 
age noise. The best commercially available microwave 
sources achieve L(Aui) = —170 dBc/Hz. The effective 
temperature T r of the pumped resonator can be calcu- 
lated as 



T r — T r o 



TT 2 

u mw -^qL(Ao;)/10 



(13) 



where T r Q is the temperature of the resonator with- 
out pumping and Z r = L/C is the characteristic 
impedance of the resonator. Now, both terms in Eq. © 



FIG. 2: Energy-level diagram of a high-frequency resonator 
and low-frequency mechanical oscillator with different possi- 
ble transitions. The first and second quantum numbers rep- 
resent the number of excitations in the resonator and me- 
chanical oscillator, respectively. The vertical arrows represent 
the environment-induced decay in the resonator, and slanted 
arrows represent driving-induced transitions. Decay in the 
mechanical oscillator is assumed to be negligibly small. 



depend on the pumping power. The first one increases 
with pumping power while the second one decreases. 
Since the highest cooling power is expected fop& T r < 1 
the effective temperature of the mechanical oscillator is 
higher than 



T, 



mO 



TT 2 

u mw _-^qL(Alj)/10 



2k B Z r 



(14) 



Thus, for microwave resonators the first term in Eq. ([9]) 
becomes important, especially for the cooling of mechani- 
cal oscillators with high resonant frequencies approaching 
the GHz range. The minimal temperature T m o is directly 
proportional to the pumping power in the limit Vo <C 1, 
which is the relevant limit in order to determine the low- 
est achievable temperature of the mechanical oscillator 
cooled by sideband cooling. 



III. QUANTUM APPROACH 

In order to achieve the quantum regime of the me- 
chanical oscillator, the temperature should be lower than 
the energy of quantum fluctuations which, together with 
Eq. (|11|) . imply the inequality 



T 



< T* < 



hui„ 



Thus the microwave resonator should be in the quan- 
tum regime as well, and the classical description is no 
longer valid. Therefore, we now turn to the analysis of 
this problem using the quantum description when the 
resonator's frequency is higher than its temperature and 
the resonator is in its ground state with high probability. 
The analysis is also valid if the resonator is substituted 
by a two-level system (qubit) as sugested in Ref. [l(J In 
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this case the cooling limit can be derived in a transparent 
manner. 

The Hamiltonian that we shall use in our analysis is 
given by 

H = LU r a\a r + LO m a\ n a m + ^coupling + #drivc, (15) 

where a \ and a r (aj^ and a m ) are, respectively, the cre- 
ation and annihilation operators of the resonator (os- 
cillator). The term H coup n ng represents the oscillator- 
resonator coupling, and the term -f/drivc represents the 
driving force. We shall assume that the last two terms 
in the Hamiltonian are small: The smallness of -ffcoupiing 
means that the energy eigenstates will, to a good approx- 
imation, be identified with well-defined excitation num- 
bers in the oscillator and resonator, while the smallness 



of -ffdrive justifies a description of the system using time- 
independent energy levels. In the following we start by 
using thermodynamics arguments to derive an expression 
for the lower limit on the achievable temperature, and we 
later use a master-equation approach to treat the specific 
example discussed in Sec. HU 

We first consider the situation depicted in Fig. [5] Each 
arrow describes a transition from a state \i,j) to an- 
other state where the meaning of the quantum 
numbers is explained in Fig. [2j We denote the rate at 
which such a transition occurs by Wuj^u/^ty In other 
words, the probability current of the transition is given 
by P\i,j)Wuj\-+\ii i ji\, where P\i,j) is the occupation prob- 
ability of the state \i, j). In the steady state, we can write 
detailed-balance equations for the occupation probabili- 
ties of the different quantum states in the form 



dPu 



dt 



{W\i+i,3-i)^\i,j)P\i+i,j-i) - W\ij)^\ i+ ij-i)P\ itj) ) + O^-ij+i^i^iVij+i) - Wiij)^\i^i tj+ i } P] itj} ) 

+ i W \i+l.j)^\i,j)P\i+l,]) - W \i,j)^\i+l,j) P \i,j)) + ( W \i-lj)->\i,j) P \i-l,3) ~ W \i,0)^\i-l-j) P \i,0)) ( 16 ) 

I 



Now we determine some relations among rates W . Let 
us start with situation when the driving force is switched 
off and the resonator is in contact with its surrounding 
environment, which is at temperature T r . Assuming that 
the environment induces transitions between states that 
are different by one photon in the resonator, the rates 
must obey the thermal-equilibrium relation 



W \i+l,j)^\i,3) 



exp 



fko r 
kuT r 



(17) 



Note that these transitions do not change the state of 
the mechanical oscillator, since without the driving force 
the oscillator and resonator are effectively decoupled 
(u) m -C uj r ). The oscillator is itself in contact with its 
environment at temperature T m , but for optimal cooling 
we assume that the insulation is good enough that we 
can completely neglect environment-induced transitions, 
i.e. we have assumed that Wuj\->\i j±i\ — ► in Fig.[2]and 
Eq. (fT6|) . We now assume that the driving force couples 
states of the form \i,j) and \i + l,j — 1) but does not 
drive any other transitions. Since the driving force is a 
classical one, the transitions it induces must have equal 
rates in both directions, i.e. 



ir. 



(18) 



The reason why there is no Boltzmann factor in Eq. (|18p 
is that these transitions are mainly induced by the clas- 
sical driving force, and any contributions to their rates 
from the thermal environment are negligible. 

Using Eqs. (fT?]) and (fTS)) . it is not difficult to verify 



that the pairs of terms in Eq. (|16[) all vanish when 

(i + j)fko r 



P \i,j) = 



z exp 



knT r 



(19) 



where Z is the partition function. This steady-state prob- 
ability distribution Pwj) can now be rewritten as 



P \i,3) 




f ihco r ] 
I k B T r ) 




exp<j 


\ jHu; r 1 

I k B T r j 




= i ex H 


f ifiw r 1 
I k B T r \ 




exp<j 


r jhw m \ 



with 



711 



(20) 



(21) 



Here Z r and Z m is partition function of resonator and 
mechanical oscillator, respectively. We therefore find 
that if the above picture about the allowed transitions 
and the relations governing their rates are valid, we can 
reach the final temperature given by Eq. (f2"Tj) . 

The above derivation suggests an intuitive picture for 
the cooling mechanism. The purpose of the driving force 
is to facilitate the transfer of excitations between the 
resonator and oscillator. Before the driving starts, the 
low-frequency oscillator has many more excitations than 
the high-frequency resonator. Once the driving starts, 
the excitation imbalance causes excitations to start flow- 
ing from the oscillator to the resonator. As the number 
of excitations in the resonator goes above the thermal- 
equilibrium value, excitations start to dissipate from the 
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FIG. 3: Same as in Fig. [2] but with the driving force inducing 
a different type of transitions. The driven transitions in this 
case remove one excitation from the oscillator state and add 
two excitations to the resonator state, or vice versa. 



Note that this situation would require driving the system 
at a frequency 2u> r — ui m . The argument can be general- 
ized to obtain any value n in the denominator of Eq. (|22j) , 
or more explicitly 



m 



n u) r 



(23) 



The question is whether such transitions can be realis- 
tically driven in a given system with a specific form of 
resonator-oscillator coupling and a given type of driving 
force. 

In order to determine the feasibility of realizing con- 
ditions where Eq. (|2"3"|) with n > 1 is the relevant lower 
limit on the achievable temperature, we now turn from 
the above general arguments to the specific situation con- 
sidered in Sec.|nl Equations ([l]),© result from a Hamil- 
tonian of the form 



resonator to the environment. A steady state is eventu- 
ally reached with the resonator in thermal equilibrium 
with the environment and both the resonator and the os- 
cillator having the same average number of excitations 
(in fact, the resonator and the oscillator will have the 
same excitation- number probability distribution). This 
picture of the cooling mechanism reveals another point 
that is generally not noted in the literature. Although T r 
is desired to be smaller than Lu m in order to avoid heating 
effects, it should not be too small, because it provides the 
mechanism by which excitations are dissipated from the 
resonator into the environment. In particular, it must 
be larger than T m , such that the dissipation of excita- 
tions is faster than the heating of the oscillator by its 
environment. 

We now consider what would happen if one were able 
to drive transitions as shown in Fig. [31 With optimal 
parameters for cooling, one would obtain the minimum 
temperature 



H = uj r a},a r + u m a* m a m + g (a r + oj) (a m + a] n ) 
+ A cos(w p t + 6) (a r +al), (24) 

where g is the oscillator-resonator coupling strength, and 
A is the amplitude of the driving force. Starting with 
a simplified version of the above Hamiltonian that does 
not contain the last two terms, we have an energy-level 
diagram similar to the one shown in Fig. [2] (without the 
induced-transition arrows). The resonator-oscillator cou- 
pling term mixes the different quantum states together in 
the eigenstates of the Hamiltonian. This mixing allows 
the driving term, which would normally affect only the 
resonator, to drive transitions that remove excitations 
from the oscillator and add excitations to the resonator, 
or vice versa. Therefore, in order to determine the transi- 
tions that can be driven by the external force, we need to 
evaluate matrix elements of the form (4>ij\(a r +al)\ipk,i) , 
where we now use the eigenstates of the Hamiltonian (i.e. 
slightly modified from the case of two uncoupled oscilla- 
tors). To first order in perturbation theory, 



\M » + {2t t 1)9 (VJ\hi - 1) - VJ+i\iJ + 1)) 

+ 7T^ (V*(i-l)0' + l)|i - 2, j + 1) - y/(i + l)(i + 2)j\i + 2, j - 1)) 



(25) 



Using the above approximation, we find that It is straightforward to see from Eq. |24)) that 

t 2 gv ^TT)7 (^K + at|^ +2lj -i) = 0. 

{iPi,j\ar + al-m+i,]-i) ~ 1 ■ (26) 
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Using numerical calculations we find that 



h3 I 



%+3,j-l) 



12g 3 y/(i + l)(i + 2)(i + 3)j 



(2u r 



0J„ 



The above results imply that the driving term can be 
used to drive transitions of the form <-> \i + 1, j — 1), 
which can be used to remove excitations from the oscil- 
lator and add them to the resonator. These transitions 
correspond to the picture shown in Fig. [2 and their res- 
onance frequency is given by lo p — Lu r — uj m . The steady- 
state effective temperature for the oscillator is given by 
Eq. (|21[) when driving these transitions, assuming that 
heating effects are avoided. By driving the system at the 
frequency uj p = 3w r — uj m , one could in principle drive 
the transitions <-* \i + 3, j — 1} and reach a lower 
minimum temperature. However, the fact that the cor- 
responding matrix element is proportional to the third 
power of the small coupling strength g suggests that this 
matrix element will be extremely small for any realis- 
tic parameters, hindering the possibility of utilizing this 
cooling mechanism. 

We now turn to the heating effects that have been ne- 
glected above. We note that the driving term in Eq. 
can also drive transitions of the form \i,j) <-> + j + 
and the relevant matrix element is given by 



«4 iv* 



i+ij+i/ 



2gy/(i + l)(j + l) 



(28) 



These undesired transitions are induced if either the driv- 
ing amplitude A or the resonator's damping rate r r 
is comparable to or larger than ui m . If either one or 
both of the above conditions are satisfied, the driving 
force must be considered within the resonance region 
of the above transition. As a result, additional excita- 
tion would be steadily pumped into the system, resulting 
in a higher temperature than what would be obtained 
from the simple picture of transition rates that we have 
presented above. The ideal parameters for cooling are 
therefore given by A <c {ijj r — Wp) and 

r r <C {uj r — uj p )\ naturally T m is desired to be much 
smaller than both T r and A 2 /T r , such that the oscillator 
heating from its contact with the environment is slower 
than the cooling it experiences as a result of the driving. 
Using a numerical simulation, we shall see shortly that 
the above heating effects can be made negligible with 
the proper choice of parameters. We should also men- 
tion here that in this section we have not considered the 
noise in the driving force, i.e. we have assumed a per- 
fect microwave source. Such noise would directly heat 
the resonator, resulting in a higher base temperature, as 
discussed in Sec. [TTJ 



In order to give a concrete example that illustrates 
the cooling dynamics, we now turn to a master-equation 
approach. The density matrix p of the system evolves in 
time according to the master equation 



dp 



+ (l + N r ) T r I a r pal - -a\.a r p ~ T,P a \ a r ) + N r T r ( alpa r - ~a r alp - -pa r a\ 



+ (1 + N m ) T m a m pa. 



1 



1 



a J m a m p - -pa] n a m I + N m T rn I a] n pa m - - a mal n p - -pa m a 



1 



1 



(29) 



where 



N,. 



phuj r fk£ 



i 



(30) 



between the effective temperatures of the oscillator and 
the resonator is approximately equal to ui m /u) r . 



and similarly for N m . The coefficients T r and T m are 
decay rates for the resonator and oscillator, respectively. 

An example illustrating the dynamics of cooling the 
mechanical oscillator by the microwave resonator is 
shown in Fig. [H The results were obtained by nu- 
merically solving Eq. (|29p using the Hamiltonian in 
Eq. (|24p . The effective temperatures of the oscilla- 
tor and resonator are obtained by calculating their re- 
spective entropies from their reduced density matrices 
(S = — Trace {p log p}) and fitting these values to the 
temperature-entropy relation for a harmonic oscillator. 
The initial heating of the resonator is a result of the trans- 
fer of excitations from the oscillator to the resonator. For 
large t, the system reaches a steady state where the ratio 



IV. CONCLUSIONS 



We have shown that both the classical and the quan- 
tum treatment give the same final result: the cooling 
factor T^/T r is limited by the ratio u) m /u) r . This lower 
limit for the cooling becomes crucial for rf and microwave 
resonators pumped by a real (noisy) microwave source 
since their effective temperature T r is usually much larger 
than the ambient temperature. We should also empha- 
size that our results apply, with minor modifications, to 
other types of coolers, e.g. a cooper-pair box. 
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FIG. 4: Effective temperatures for the resonator (solid line) 
and the oscillator (dashed line), relative to the ambient tem- 
perature, as a function of time for the parameters uj m juj r — 
0.1, g/ui r = 0.005, k B T m .rlhu] r = 0.5, r r /(w r /27r) = 0.05, 
and r m = 0. The driving field, with amplitude A/u r = 0.05 



and frequency lu p — uj t 



is turned on at t = 0. 
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